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^ I Abstract 

^ ■ In this paper we obtain the Fourier Transform of a continuous gravitational wave. 

O We have analysed the data set for (i) one year observation time and (ii) arbitrary obser- 

' vation time, for arbitrary location of detector and source taking into account the effects 

arising due to rotational as well as orbital motion of the earth. As an application of the 
VHy transform we considered spin down and N-component signal analysis. 



1 Introduction 



Detection of gravitational waves (GW) mainly depends on the efficient Fourier analysis of 
the data of the output of the detector such as LIGO/VIRGO. In an earlier paper we presented 
Fourier analysis of one day observation data set of the response of a laser interferometer 
(Srivastava and Sahay, 2002 a). Hereafter the reference to this paper and its results are 
noted by the prefix I. We have seen that the amplitude and the frequency modulations 
result into a large number of side bands about the signal frequency f . Consequently, 
the maximum power lies in the frequency f a + 2f rot with amplitude reduction by 74% to 
what one would have expected due to increased data interval. Here f rot is the rotational 
frequency of the earth. Hence, for GW detection it is desirable to obtain FT for longer 
observation data. In the next section we investigate one year observation data. Incidently 
there exists correspondences and identifications of the present analysis to the analysis of 
paper-I. To facilitate analogous modifications we introduce corresponding quantities with 
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tilde viz., C and V in place of C and V. We obtain in section 3 generalisations of the 
results for arbitrary observation time. As an application of the results obtained we consider 
spin down and N-component signal analysis in sections 4 and 5. Last section contains the 
discussion and conclusions of the paper. 



2 Fourier transform for one year integration 

The phase of the GW signal of frequency f a at time t is given via (1-25) i.e. 



S(f) 
where 
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where R e , R se , w rot and w orb represent respectively the earth's radius, average distance of 
earth's centre from the origin of Solar System Barycentre (SSB) frame, the rotational and 
the orbital angular velocity of the earth. 9 and denote the celestial colatitude and celestial 
longitude of the source, e and c represent the obliquity of the ecliptic and the velocity of 
light, a is the colatitude of the detector. Here t represents the time in seconds elapsed from 
the instant the sun is at the Vernal Equinox and (3 Q is the local sidereal time at that instant, 
expressed in radians. 



The Fourier transform (FT) for one year observation time, T obs is given as 



Kf) 



aT 



cos[$(t)]e- i2 ^t ; 



y Jo 

a = a~ l = w rot /w or b ; T = one sidereal day ; 
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T obs = aT ~ 3.14 x 10 7 sec. (6) 
This splits, similar to Eqs. (1-33-35), into two terms as 

Kf)].=h-+h + \ (7) 



j_ = 1 f 2n e i[£i>-+Zc™(Z-4>)+Mcos(aZ-5)-n-Q] d F ^ 

2w or b Jo 



j _ 1 f 2w e -i[£v + +Zcos(£-4>)+Arcos(a£-5)-Tl-Q] c i£ 

2w orb Jo 



(9) 
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V T = —L 5 f = ^orb = W orb t (10) 

Jorb 

where / or b is the orbital frequency of the earth. Hereafter, we neglect I p+ as it oscillates 
rapidly and contributes very little to h(f) and write v in place of P_. A careful compar- 
ison of Eq. (§) with (1-33) i.e. 
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obtained for one day observation data reveals that the integrand of the equations are identi- 
cal with following identifications and correspondences. 

5 *-* 4> , 

Z «-> AT, \ (12) 
a <-> a. 

We may, therefore, employ the results obtained there by introducing obvious corresponding 
quantities i.e. B, C, V in place of B, C, V leaving A unchanged. Hence 

_ — k=oo m=oo 

~ h (f)] v ^^— E E e^B[C-W); (13) 

y ^UJorb k=-oo m=-oo 

A = tegte-K-Q, 

g _ J k (M)Jm(Z) 

v 2 — (ak+m) 2 ' 

C = sin 2z/7r cos(2a/c7r — kS — m</>) — afc + rn | cos 2utt sin(2a/c7r — k5 — rruf)) + sin(k5 + m0)} 
X> = cos 2vi\ cos(2a/c7r — k5 — m<f)) + ka + m sin 2P7r sin(2a/c7r — kS — m0) — cos(M + m</>) 

where J stands for the Bessel function of first kind. Now the FT of the two polarisation 
states can be written as 
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The FT obtained contains double series of the Bessel functions of the orders k and m 
ranging from — oo to oo. It is well known that the Bessel functions decrease rapidly as the 
order exceeds the argument. Hence possible range of k and m over which the summation 
of the series is to be considered depends on the arguments of Bessel functions i.e Z and J\f. 
Referring to Eq. O) it is found that 
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The FT of a FM signal for 
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is shown in Fig. The spectrum resolution is kept equal to 1/T Q i.e. 3.17 x 1CT 8 Hz. 
We have convinced ourselves by plotting the FT at higher resolutions that the resolution 
of 1/T is sufficient to represent the relevant peaks. We notice that the drop in amplitude 
at the source frequency is about 98% which may be attributed to the presence of a large 
number of side bands. 

The complete response R(t) of the detector and its FT can be obtained similar to the 
one presented in section 4 of the paper-I. Taking the results straight away from (1-43, 44, 
45) one may get 
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Figure @ shows the power spectrum of the complete response of the Doppler mod- 
ulated signal. We have kept here all the parameters same as in FM. In this caee also we 
observe that maximum power is associated with f Q + 2f rot and the least lies in f Q — f rot . 

3 Fourier transform for arbitrary observation time 

It is important to obtain the FT for arbitrary observation time. The results obtained will 
be employed to outline how spin down of a pulsar due to the gravitational radiation back 
reaction or due to some other mechanism can be taken into account. 
The FT for a data of observation time T is given via 

h(f) = [ T ° cosmt)}e- l2nft dt (20) 
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which may be split 
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[refer to Eqs. I (33 - 35)]. 

As I v+ contributes very little to h(f) we drop I u+ and write v in place of z/_. Using the 
identity 
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After performing the integration and proceeding in a straight-forward manner we have 
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The FT of the two polarisation states of the wave can now be written as 
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Now it is simple matter to obtain the FT of complete response. One gets 
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The FT of FM signal of a detector for a data of 120 days for 

f = 25 Hz, h 0+ = h 0x =1 j 
« = 7r/6, /3 d = 7t/3, 7 = 27r/3, 
6 = tt/9, = tt/4, ^ = tt/4. 



(3D 



(32) 



is plotted in Fig. @ with a resolution of 1/2T Q ps 9.67 x 1CT 8 Hz. The Power spectrum of 
the complete response is plotted in Fig. (Q). In this case also the most of the power of the 
signal lies in the frequency f a + 2f rot and the least with f a — f rot . 
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4 Spin down 



Pulsars loose their rotational energy by processes like electro-magnetic breaking, emission 
of particles and emission of GW. Thus, the rotational frequency is not completely stable, but 
varies over a time scale which is of the order of the age of the pulsar. Typically, younger 
pulsars have the largest spin down rates. Current observations suggest that spin down is 
primarily due to electro-magnetic breaking (Manchester, 1992 and Kulkarni, 1992). Over 
the entire observing time, T Q the frequency drift would be small but it may be taken into 
account for better sensitivity. To account this aspect we consider the evaluation of FT in a 
sequence of time windows by splitting the interval — T Q in M equal parts, each of interval 
At (T Q = M At) such that the signal over a window may be treated as monochromatic. The 
strategy is to evaluate the FT over the window and finally to add the result. This process has 
been suggested by Brady and Creighton (2000) and Schutz (1998) in numerical computing 
and has been termed as stacking and tracking. For any such window the initial and final 
times may be taken respectively as t Q + n A t and t Q + (n + 1) A t where t Q represents the 
starting of the data set and < n < M — 1. The window under consideration is the n th 
window. Let 

r t +(n+l)At 
I = / h(i)e- t2nft dt 
Jt +nAt 

= f At hit + t + n A t)e- t2w{t+to+nAt)f dt ; (33) 
Jo 

t = t + t + nAt (34) 
Hence the FT for the data of the time-window under consideration is given via 

'h(f)] = [ At CO s[$(t + t + n At)]e- l2 ^ t+t ° +nAt ^ f dt (35) 
L is Jo 

Taking the initial time of the data set 

to = (36) 
and proceeding as in the previous section we obtain 
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After integration we get 
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A s = ^^ + 2nnAt(f -f)-n-Q 

g = J k {Z)J m {N) 

v 2 — (afc+m) 2 ' 

C s = sin(z/r) cos(akr + mr — k\ — mC) — afc + m {cos(z/r) sin(a£;r + mr — k\ — m() 
+ sin(fcA + m£)} , 

P s = cos(z/r) cos(akr + mr — k\ — mQ + afc+m sin(z/r) sin(a£;r + mr — k\ — mQ 
— cos{k\ + to£) , 

A = <p — aw , £ = 5 — nr , 

r = uwAt, n = 0,1,2,3, ,M-1. 

The FT of the complete response would now be given via 
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5 N-component signal 

The FT in Eqs. ( p7| ) and ( |3T| ) are valid for a pulsar which emits GW signal at single fre- 
quency /„. But there are known physical mechanisms which generate GW signals consist- 
ing of many components. An axially symmetric pulsar undergoing free precession, emits 
quadrupole GW at two frequencies, one equal to the sum of the spin frequency and the 
precession frequency, and the other twice of it (Zimmermann and Szedenits, 1979). The 
quadrupole GW from a triaxial ellipsoid rotating about one of its principal axes consists of 
one component only (Thorne, 1987). In this case the signal has frequency twice the spin 
frequency of the star. In general, if a star is non-axisymmetric and precesses, the GW signal 
consists of more than two components. Recently, new mechanisms e.g. r-mode instabil- 
ity of spinning neutron stars (Anderson, 1998; Lindblom et al., 1998; Owen, et al., 1998) 
and temperature asymmetry in the interior of the neutron star with miss-aligned spin axis 
(Bildsten, 1988) have been discussed in the literature. 

In view of the above discussion CGW signal may consist of frequencies which are 
multiple of some basic frequencies. An analysis of the GW data of N-component of signal 
has been made recently by Jaranowski and Krolak (2000). We do not continue this analysis 
as the requisite formalism is analogous to what we have presented in paper-I and the earlier 
sections of this paper 
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6 Discussion and conclusions 



The analysis and results obtained in paper-I regarding FT of the response of a Laser In- 
terferometer have been generalised in this paper. In this context following points must be 
noted. 

1. For 120 days observation time data the resolution provided by FFT, which is equal 
to 1/T , is sufficient to represent the structure of side bands. 

2. In every case discussed, it turned out that the maximum power lies in the frequency 
f + 2f rot . However, this is not established conclusively if this result is generic. 

3. Throughout our analysis in paper-I and this paper we have employed following con- 
ditions. 

(i) The phase of the wave is zero at t = 0. 

(ii) The observation time of the data set is from t = to t = T Q . 

4. By trivial modifications we could get the general results with non-zero initial time 
and phase. 

It is pointed out that in choosing the parameters for the source and the detectors we have 
been constrained by the computer memory and efficiency. Accordingly the choice of the 
parameters has been a matter of convenience. The application of the results obtained in the 
present and the earlier paper to the problem of all sky search for continuous gravitational 
wave sources will be made in our forthcoming article (Srivastava and Sahay, 2002 b). 
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Figure 2: Power spectrum of the complete response of a modulated signal of frequency, 
f = 50 Hz from a source located at (7r/18, 0) with a resolution of 3.17 x 1CT 8 . 
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Figure 3: FT of a FM signal of frequency, f = 25 Hz from a source located at (n/9, 7r/4) 
with a resolution of 9.67 x 10~ 8 . 
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Figure 4: Power spectrum of the complete response of a Doppler modulated signal of 
frequency, f = 25 Hz from a source located at (n/9, 7r/4) with a resolution of 9.67 x 1CT 8 . 
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